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Abstract. In creating his gravitational field equations Einstein unjustifiedly assumed that inertial mass, 
and its energy equivalent, is a source of gravity. Denying this assumption allows modifying the field equations 
to a form in which a positive cosmological constant appears as a uniform density of gravitationally repulsive 
matter. This repulsive matter is identified as the back sides of the 'drainholes' (called by some 'traversable 
wormholes') introduced by the author in 1973, which attract on the high, front sides and repel more strongly 
on the low, back sides. The field equations with a scalar field added produce cosmological models that 
'bounce' off a positive minimum of the scale factor and accelerate throughout history. The 'dark drainholes' 
that radiate nothing visible are hypothesized to constitute the 'dark matter' inferred from observation, 
their excess of negative active mass over positive active mass driving the accelerating expansion. For a 
universe with spatial curvature zero, and the ratio of scale factor now to scale factor at bounce equal to the 
Hubble radius over the Planck length, the model gives an elapsed time since the bounce of two trillion years. 
The solutions for negative spatial curvature exhibit early stage inflation of great magnitude in short times. 
Cosmic voids, filaments, and walls are attributed to separation of the back sides of the drainholes from the 
front, driven by their mutual attractive-repulsive interactions. 
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Cosmologists are perplexed by the discovery that the expansion of the universe, long thought to be 
slowing down, is in fact speeding up. This acceleration seems to require, in addition to the mysterious, 
unseen 'dark matter' invoked to explain the assembling of visible matter into galaxies and galactic clusters 
and superclusters, an even more mysterious 'dark energy' that acts in a gravitationally repulsive manner to 
cause the acceleration. The initial attempt to model this dark energy by Einstein's cosmological constant 
A, attributing its source to a negative pressure created in the vacuum by virtual particle pairs, runs up 
against a discrepancy of at least fifty-five orders of magnitude between the tiny A required to explain the 
acceleration and the large A predicted by the quantum mechanics of these virtual pairs. A way out of this 
conundrum is available, but it requires the correcting of an erroneous assumption Einstein made in the early 
days of the general theory of relativity, an assumption that has stood virtually unchallenged down through 
the years. The assumption in question is not the introduction in 1917 of A into the field equations of gravity, 
self-described by Einstein as a mistake. It had in fact already appeared in his 1916 paper Die Grundlage der 
allgemeinen Relativitatstheorie pQ setting out the fundamentals of the general theory. 



Selected for honorable mention in the Gravity Research Foundation 2006 Awards for Essays on Gravitation. 
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In its most transparent form the assumption is that inertial mass, and concomitantly its energy equiva- 
lent, is a source of gravity and must therefore be coupled to the gravitational potential in the field equations 
of the general theory. Einstein arrived at this assumption in the 1916 paper while seeking a tensorial equation 
to correspond to the Poisson equation V 2 </> = Aitk/j,, where \i denotes the "density of matter" . Drawing on the 
special theory's identification of "inert mass" with "energy, which finds its complete mathematical expression 
in . . . the energy-tensor" , he concluded that "we must introduce a corresponding energy-tensor of matter 
T°" . Further describing this energy-tensor as "corresponding to the density /j in Poisson's equation" , he 
wrote down his now hallowed field equations R a p — \Rg a [3 — ^TrT a p. The unjustified step in this argument 
is the confusing of 'gravitating mass', which is the sole contributor to the "density of matter" in Poisson's 
equation, with "inert mass", which is indeed equivalent to energy in the proportion m = E/c 2 . That all 
bodies respond alike to a gravitational field establishes the equivalence of 'inertial' (inert) mass with 'passive' 
(gravitated) mass, but there is no corresponding link between passive and 'active' (gravitating) mass, thus 
no link between inertial mass and active mass. In Newton's theory of gravity there is such a link, but it 
depends on his law of action and reaction, which for gravitating bodies would require instantaneous action 
at a distance, something that in a relativistic field theory such as Einstein's does not exist. It likely is 
this inference from Newton's theory that caused Einstein to treat "inert mass" and "density of matter" as 
equivalent. 

If, contrary to Einstein's assumption, inertial mass and its energy equivalent is not a source of gravity, 
then in particular the kinetic energy in the form of the pressure p in a continuous distribution of gravitating 
matter must not contribute to gravity, so Einstein's choice T Q/3 = pu a u 13 + (p/c 2 )(u a u^ — g a/3 ) must be 
modified. The most elegant way to effect a modification is to not think about T a p, but instead derive field 
equations as the Euler-Lagrange equations of the action integral 

j{R-^){-g)hd , X, (1) 

where /x is the active gravitational mass density. (This integral is the most straightforward relativistic analog 
of the action integral that yields the Poisson equation for the newtonian gravitational potential.) Variation 
of the metric produces the modified field equations 

R a0 -^Rg a f3 = -^^g a0 , (2) 

which makes T a p — —\^g a p- Equivalent to Eq. (2) is R a p — ^rf-ig a p, the 00 component of which reduces 
in the slow motion, weak field approximation precisely to the Poisson equation. 

Incorporation of other putative determinants of the geometry of space-time, such as scalar fields and 
electromagnetic fields, can be accomplished in the usual way by adding terms to the action integrand. In 
particular, a cosmological constant term can be added, changing the integrand to R — ^§rp + 2A and the 
field equations to 

Rap-\Rg*p = -^{v + n)g a p, (3) 

where ^rp, = —A. Seen in this light a positive A is simply a (mis)representation of a negative active 
mass density p of a continuous distribution of gravitationally repulsive matter, an excess of which over the 
positive active mass density /j of attractive matter could drive an accelerating cosmic expansion. Here one 
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sees a glimmer of a solution to the 'Cosmological Constant Problem'. The question is: Why should such 
gravitationally repulsive matter exist and where should we look for it? 

In 1973 I described in considerable detail a model of a gravitating particle alternative to the Schwarzschild 
vacuum solution of Einstein's field equations. This space-time manifold, which I termed a 'drainholc', has 
subsequently come to be recognized as an early (perhaps the earliest) example of what is now called by 
some a 'traversable wormhole' [2j |4]- The metric is a static, spherically symmetric solution of the field 
equations R a p — ^Rg a p = T a p :— —2(<f>. a <f>.p — ^</>' 7 <^. 7 g a p) and := </> 7 :7 = arising from the action 
integrand R + 20 7 7 . (N.B. R a p and R here are the negatives of those in [2].) It has the proper-time form 
(in units in which c = 1) 

dr 2 = [1 - f(p)} dT 2 [1 - f(p)}- 1 dp 2 r 2 (p) dn 2 

(4) 

= dt 2 - [dp - f(p) dt} 2 - r 2 (p) dn 2 , 
where t = T-J f(p)[l - P(p)]- 1 dp, 

r(p) = ^J(p - m) 2 + a 2 e (™/«)«(p) and 1 _ = e -{2m/n) a ( P ) ^ ^ 



a(p) = - 
a 



7T 

tan 

2 



(G) 



and a := y/n 2 — to 2 , to and n being parameters satisfying < m < n. (The coordinate p used here is the p 
of [2] shifted upward by m.) The shapes and asymptotics of r and f 2 are shown in Fig. 1. Not shown, but 
verifiable, is that f 2 {p) ~ 2m/ p as p — > oo. 




/ 2 (P) f 



— 2mn/ a 




2m 



Fig. 1. Graphs of r(p) and f 2 (p) for typical values of the parameters to and n. 

The choke point of the drainhole throat is the 2-sphere at p — 2m, of superficial radius r(2m) (i.e., of 
surface area 47rr 2 (2m)), which increases monotonically from n to ne as m increases from to n. Thus the size 
of the throat is determined almost exclusively by n, independently of m. Although the scalar field <p has a 
nonkinctic 'energy' density that contributes to the space-time curvature through T a p, this energy has little to 
do with the strength of gravity (as determined by m) , rather is associated with the negative spatial curvatures 
found in the open throat, the negativity of which mandates the minus sign at the front of T a p. Because 
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r(p) > n > and f 2 (p) < 1, the space-time manifold is geodesically complete and has no one-way event 
horizon, the throat being therefore traversable by test particles and light in both directions. The manifold 
is asymptotic as p — > oo to a Schwarzschild manifold with (active gravitational) mass parameter m. The 
flowing 'ether' (a figurative term for a cloud of incrtial observers free-falling geodesically from rest at p = oo) 
has radial velocity f(p) (taken as the negative square root of f 2 (p)) and radial acceleration (f 2 /2)'(p), which 
computes to —m/r 2 (p) and therefore is strongest at p = 2m. Because the radial acceleration is everywhere 
less than 0, the drainhole attracts test particles on the high, front side, where p > 2m, and repels them on 
the low, back side, where p < 2m. Moreover, the manifold is asymptotic as p — > — oo to a Schwarzschild 
manifold with mass parameter m = -me™'", so the drainhole repels test particles more strongly on the low 
side than it attracts them on the high side, in the ratio —m/m = e m7r /Vn 2 -m 2 _ rpj^ drainhole [ s a kind f 
natural accelerator of the 'gravitational ether', drawing it in on the high side and expelling it more forcefully 
on the low side, much as a leaf blower does air. 

The 1973 paper has in it the following sentence: "A speculative extrapolation from the asymmetry 
between m and m is that the universe expands because it contains more negative mass than positive, each 
half-particle of positive mass m being slightly overbalanced by a half-particle of negative mass m such that 
— m > m." To bring this idea to bear on the cosmological conumdrum, let us study solutions of field 
equations that incorporate a positive mass density p, a negative mass density p such that —p > p, and a 
scalar field <p as above. The action integrand R — ^f{p + p) + 2</r 7 </>. 7 combines these elements and yields 
the field equations 

Rap - \Rg a f) = T af} := -^£(m + P) g a p - 2{<f>. a (p.0 - l<j>-~<<j>. 7 g a i3) (7) 

and D(j) := </r 7 :7 = 0. For a Robertson- Walker metric c 2 dt 2 — R 2 (t)ds 2 and a scalar field = (3(t) these 
reduce to 

n R 2 /c 2 + k 4™ , p 2 



^R 



R 2 /c 2 



R 2 



Ann 



■P) 



ff 2 



(9) 



and 



□0 = ? + 3/3- J = 0, (10) 

where k = 1, 0, or —1, the uniform curvature of the spatial metric ds 2 . In addition there is, corresponding 
to the identity Tj 3 .^ = 0, the equation ^£-d(p + p) = —2(\3<j))d<j) = 0, which implies that the 'accelerant' 
A, defined by A := —^r{p + p), is a constant, positive under the assumption that — p, > p. Equation (10) 
integrates to $ 2 R 6 = Be 2 , where B also is a positive constant if $ ^ 0. Equations (8) and (9) then arc 
together equivalent to 

1 R 2 _ 4ttk k $ 2 _ A k B _AR 6 -3kR i -B 

^R? ~~ ~lk?^ + ^> ~ R? ~ 3^2 - 3 ~ & ~ 3i?e - 3~R6 (U) 



EXPLAIN COSMIC ACCELERATION? FIRST, CORRECT EINSTEIN 



5 



and 



1 R 4ttk, _ n 2/3 2 A 2B AR 6 + 2B 



'R = -^ + ^ + ^ = ^ + 3Rt = ^R^- (12) 

Four implications of these equations are immediate. First, the scale factor R has a positive minimum 
value i? m in (namely, the only positive root of the polynomial AR 6 — 3kR 4 — B), which rules out a 'big 
bang' singularity, putting in its stead a 'bounce' off a state of maximum compression at time t = 0, when 
R(t) = i?min- Also, R(t) — > oo as t — > ±oo. Second, R is always positive, so the universal expansion is 
accelerating at all times after the bounce, and the universal contraction is decelerating at all times before 
the bounce. Third, the 'Hubble parameter' H (:= R/R) behaves asymptotically as follows: 

1 2 A 3kR 4 + B A [from below if k > ol , N 

-~H = 1 > - I ~ } as R -> oo. 13 

c 2 3 3R 6 3 [from above if k < Oj V ; 

Fourth, the 'acceleration parameter' Q (:= (R/R)/(R/R) 2 ) behaves this way: 

^ , kR 4 + B _ [from above if k > ol 

Q = H 5 > 1 < ~ > as i? -> oo. (14) 

F 2 i? 6 [from below if k < Oj V ' 

Lacking singularities and horizons, and being geodesically complete, mathematical drainholes are more 
pleasing to the aesthetic sense than are mathematical blackholes. Because in principle they are able to re- 
produce all the externally discernible aspects of physical blackholes that mathematical blackholes reproduce, 
they are at least as satisfactory as the latter for modeling centers of gravitational attraction. In this role 
they are more aptly called 'darkholes', inasmuch as they can capture photons that venture too close, but, 
unlike blackholes, must eventually release them, either back to the attractive high side whence they came or 
down the throat and out into the repulsive low side. Thus one can imagine that at galactic centers will be 
found not supermassive blackholes, but supermassive darkholes instead. This, however, is not the end of the 
story. A central tenet of the general theory of relativity is that every elementary object that 'has gravity' is 
a manifestation of a local departure of the geometry of space-time from flatness. If such an object has other 
properties ascribed to it by quantum mechanics or quantum field theory, these must be additional to the 
underlying geometrical structure. I therefore propose the hypothesis that every such elementary gravitating 
object is at its core an actual physical drainhole/darkhole (a 'dark drainhole' I will call it) — these objects 
to include not only elementary constituents of visible matter such as protons and neutrons, or, more likely, 
quarks, but also the unseen particles of dark matter whose existence is at present only inferential. It then 
becomes a question of to what extent this hypothesis, coupled with the cosmological model described above, 
fits the current state of observational cosmology. 

The significant parameters of the model are k, A, B, H, Q, to (the present epoch), and i? m i n (= R(0)) 
and R(t ) (or perhaps only R(to)/R m in)- Although the current suspicion that space is perfectly flat (k = 0) 
is perhaps not applicable in the context of this model, let us proceed for the moment on the presumption 
that it is. It then becomes straightforward to integrate Eqs. (11) and (12), the result being 

R 3 (t) = R^coshiVzAct), (15) 
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where i? m i n = (B/A) 1 / 6 , from which follow 



H(t) = cJ^ tanh(V3A c t) = (sgn t) c 
V o 



R(t) 



(16) 



Q{t) 1 + sinh 2 ^/3lci) 1 + [i?(t)/i? min ] 6 - 1 ' 



and 



c 2 ,4 = ff 2 (i)[Q(i) + 2] = 3H 2 (t) 



1 



(18) 



[i?(t)ARmin] 6 -l_ 

Of the present values of these parameters the only one that is reasonably well determined by observations 
is H(t ), which currently is estimated to be about 72 (km/sec)/Mpc. After H(t ) is input the others are 
determined by the ratio R{t )/R m \ n . Knowing neither the numerator, the denominator, nor the ratio itself, 
but suspecting that the ratio is quite large, the best one can do is make guesses and calculate the results. 
In this spirit let us go to extreme limits and take R(to) cm to be the Hubble radius c/H(to) (the 'radius of 
the observable universe') and R m { n cm to be the Planck length. Then R(t )/R min = 1.28 x 10 28 cm/1.62 x 
1CT 33 cm = 7.93 x 10 60 , which makes Q(t ) = 1 + 1(T 365 , c 2 A = 1.63 x lCT 35 /scc 2 = 1.62 x l(r 20 /yr 2 , and 
t = 1.91 x 10 12 years. This value for to encompasses 140 of the 13.6 x 10 9 years predicted to have elapsed 
since the 'big bang' by the 'standard' (or 'concordance') model based on the Friedmann-Robertson- Walker 
equations, an interval which in the present instance would allow approximately only a doubling from i? m i n 
to R(t). Other guesses are left to the reader, but with the advice that so long as R(to)/Rmin ^> 1 the 
parameters other than to will differ little from those values calculated above. 

Consider now the problem of estimating the ratio —p/l-i. From A = — ^r(n + p) follows — p/fj, = 
1 + c 2 A/4ttkh. If we assume, for simplicity's sake, that the dark drainholes whose active mass density is /x 
all have at each epoch the same mass and size parameters m and n, then — p/fj, = —fh/m = e m7r / a , so that 
m/\/n 2 — m 2 = m/a = ln(— Ji/ ^)/-k = ln(l + c 2 A/4ttkh)/tt. This implies that 

n 2 = [ln(l+c 2 A/47r^)] 2 

V n ) tt 2 + [ln(l + c 2 A/4ttkh)} 2 ' 1 * 

Because it is only the combination — ^pr(/U + p.) (the 'accelerant' A) that remains constant, it is possi- 
ble for the density /i to change over time in some arbitrary fashion if p, changes to compensate. Indeed 
/i (and p) might not change at all, which would indicate a 'steady-state' universe with continuous cre- 
ation of dark drainholes sustaining the expansion. In that case \x would be at all epochs the density 
at the present epoch, which according to current best estimates is the critical density /j, c of the FRW 
standard model, namely, /x c = 3H 2 (t )/8irK = 9.7 x lCT 30 g/cm 3 . From Eq. (18) one has c 2 A/Attk^ = 
ZH 2 {to) [l + l/([R(to)/R m in] 6 - 1)] /4ttk^ c = 2 [l + 1(T 366 ] = 2.0. Equation (19) then yields m/n = 0.33. 
If n is the Planck length, then m = 0.33 (1.6 x 10 -33 ) cm, = 7.2 x 10 6 grams (= 0.33 Planck mass) in 
c = k = 1 units. This makes the dark drainhole particles gravitate (not 'weigh') much more than protons 
and neutrons. To maintain the density /x c these particles would have to be created at a rate that would keep 
on average about one in every 10 9 cubic kilometers, which would keep them on average about one thousand 
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kilometers apart. To decrease the mass m and thereby increase the number density would require taking n 
smaller than the Planck length. 

Now consider the more orthodox supposition that the active gravitational mass content of the uni- 
verse is unchanging, so that \x decreases in inverse proportion to the cube of the scale factor R: thus 
/i = [Rmin/R(t)} 3 = n\t=to [R(t )/R(t)} 3 = u c [R(t Q ) I R{t)} 3 , where /x = u c [R(h) / R mia } 3 = 4.9 x 

3 



10 153 g/cm J , the density at the time of maximum compression. Equation (19) now reads 

[ln(l + (c 2 A/At:ku c ) [R(t) / R(t )] 3 )} 2 



/m 
V n 



7T 2 + [ln(l + (c 2 A/Attku c ) [R(t) / R(t )] 3 )] ' 



(20) 



The ratio of m to n increases monotonically as t goes from to oo. When t = 0, m/n = 1.3 x 1CP 183 . When 
t = to, m/n = 0.33 as in the steady-state case. As t — > oo, m/n — > 1 (the flow of the 'gravitational ether' 
through the drainholes grows asymptotically to the maximum rate that the drainholes can accommodate). 
In contrast to the steady-state version, which drives the accelerating expansion by continually producing 
new drainholes of fixed size and mass, this version drives it by continuously increasing the masses of a fixed 
population of equal-sized drainholes. The same effect could result from a mixture of the two. Also, in neither 
case is it written in stone that the sizes must be uniform — only the ratio of m to n is determinate. In 
particular, n for particle constituents of atomic nuclei should perhaps be of the order of 10 -52 cm to account 
for the implication of newtonian gravitational theory that their active masses must bear some approximate 
numerical proportion to their inertial rest masses. 
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Fig. 2. Graphs of H(R) and lnR(t), showing early stage inflation for k 



-1. 



When k = or 1 there is no early stage of extraordinary inflation, only a steadily increasing Hubble 
parameter H(t). The situation is quite different for k = —1 (strictly k = — 1/cm 2 ), as the graphs in 
Fig. 2 demonstrate. By sending B toward you can get as much inflation as you might want — and the 
more you get, the earlier you get it. One can show that i? m ; n ~ {/B/{— 3k) as B — > 0. If, to illustrate, 
A = 9.09 x 10" 57 /cm 2 and we take for R mia cm the Planck length 1.62 x 10~ 33 cm, then B w 1.94 x 
10- 131 /cm 2 , which gives H(R) the peak value H max = c^A/'i + 2 (-fc) 3 / 2 /(3 VB) = (5 x 10 m )H(t ) = 
3.6 x 10 62 (km/sec)/Mpc = 1.17 x 10 43 /sec, occurring at R = {/B/{-k) = 2.10 x 10~ 33 w 1.32 R min . 



EXPLAIN COSMIC ACCELERATION? FIRST, CORRECT EINSTEIN 



8 



Numerical solution of Eq. (12) shows that at t = 6.74 x 1CT 14 seconds the scale factor has inflated to 
R(t) = 2.02 x 10~ 3 , for a ratio R(t)/R m - m = 1.25 x 10 30 , which works out to 100 doublings. The acceleration 
Q(R) is infinite at i? m i n , is equal to 1 when H(R) = H max , bottoms out with the value Q m i n = 9.28 x 10~ 82 
at R — 4.52 x 10~ 13 , when t = 1.51 x 10 -23 seconds, then returns slowly to 1 as R —> oo (not until 
R = 1.80 x 10 29 does Q(R) reach 0.99, at which time t = 1.81 x 10 18 seconds = 5.74 x 10 10 years). Moving B 
closer to increases H max and drives i? m ; n and Q m ; n toward 0, and the 100-doublings time toward seconds. 

A final question: If dark matter, which is understood to be distributed unevenly in the universe, consists 
of the gravitationally attractive front sides of dark drainholes, where do the repulsive back sides reside? 
When some local concentration of spatial curvature (a 'quantum fluctuation', say) develops into such a 
topological hole, the entrance and the exit, if close together in the ambient space, will drift apart as the exit 
repels the entrance more strongly that the entrance attracts the exit. Apply this to a multitude of particles 
and you will likely see the front side entrances being brought together by both their mutual attractions and 
the excess repulsion from the back side exits. The exits, on the other hand, will tend to spread themselves 
more or less uniformly over regions from which they have expelled the entrances. Herein lies a mechanism 
for producing the voids, filaments, and walls of the cosmos. What is more, the filaments and walls should 
be more compacted than they would be if formed by gravitational attraction alone, for the repulsive matter 
in the voids would increase the compaction by pushing in on the clumps of attractive matter from many 
directions with a nonkinctic, positive pressure produced by repulsive gravity, not to be confused with the 
negative pseudo-pressure conjectured in the confines of Einstein's assumption to be a source of repulsive 
gravity. 

It is tempting to speculate that Einstein might have taken some of the steps described above had he 
recognized that his equating of active gravitational mass with inertial mass was unjustified. If there is a 
moral here, it would be this one, which harks back to the discovery of noneuclidean geometry that made 
possible the general theory of relativity: Examine assumptions diligently, and when you find one you can't 
justify, assume the opposite and see where it takes you. 
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